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SUMMARY 


When a low arch (a thin curved beam of small curvature) is subjected 
to a lateral loading acting toward the center of curvature, the axial 
thrust induced by the bending of the arch may cause the arch to buckle 
so that the curvature becomes suddenly reversed. The critical lateral 
loading depends on the dimensions and rigidity of the arch, the elasticity 
of the end fixation, the type of load distribution, and the initial 
curvature of the arch. A general solution of the problem is given in 
this paper, using the classical buckling criterion which is based on the 
stability with respect to infinitesimal displacements about the equilib- 
rium positions. 

For a sinusoidal arch under sinusoidal loading, the critical load 
can be expressed exactly as a simple function of the beam dimension 
parameters. For other arch shapes and load distributions, approximate 
values of the critical load can be obtained by summing a few terms of 
a rapidly converging Fourier series. The effects of initial end thrust 
and axial and lateral elastic support are discussed. 

The buckling load based on the energy criterion of Karman and Tsien 
is also calculated. The results for both the classical and the energy 
criteria are compared with experiments made on a series of centrally 
loaded, pin-ended arches. For larger values of a dimensionless param- 
eter which is proportional to the ratio of the arch rise to the 

arch thickness, the experimental critical buckling loads agreed quite 
well with the classical criterion, but, for smaller values of the 

experimental critical loads were appreciably below those calculated from 
the classical criterion, although they were always above those obtained 
from the energy criterion. 


INTRODUCTION 


An arch subjected to lateral loads may become elastically unstable. 
Generally speaking, there are two possibilities of buckling: 
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(1) If the rise of the arch (a in fig. 1) is of the same order as 
the span of the arch, then it is possible for the arch to buckle at the 
critical pressure in the mode indicated by the dashed curve in figure 1. 
Buckling of this type can be safely assumed to be ’’inextensional," as 
suggested by Lord Rayleigh, and, as such, has been discussed by 

E. Hurlbrink, E. Chwalla, R. Mayer, E. Gaber, E. L. Nicolai, and 
S. Timoshenko. (See Timoshenko’s book, reference 1, for references to 
original papers.) In all these studies, circular arches under uniformly 
distributed lateral loading are assumed, with various types of end 
fixations. 

(2) If the rise a of the arch is much smaller than the span L, 

(fig. 2), then the induced axial thrust plays an important role in the 
elastic stability. The beam may become unstable and suddenly reverse 
its curvature, jumping, for example, from the solid-line position in 
figure 2 to the dashed-line position. 

It is the object of the present paper to treat arches of small rise; 
therefore, the buckling deformation will be "extensional” rather than 
"inextensional.” It will be shown that the variation in the initial 
curvature of the beam has a very important effect on the critical load. 
Furthermore, with a view to possible applications to thin-wing design 
problems, beams acted on by initial thrust and those with elastic sup- 
ports will be discussed. 

The same problem has been treated before by Biezeno (reference 2), 
Marguerre (references 3 and 4), Timoshenko (reference 1), and Friedrichs 
(reference 5)*^ Biezeno and Timoshenko derived the fundamental dif- 
ferential equation in the same manner as this paper, while Marguerre 
and Friedrichs derived their equations by variational principles. The 
resulting equations are the same. Biezeno treated a circular arch under 
a concentrated load at the center and Marguerre and Friedrichs, a cir- 
cular arch under uniformly distributed pressure; all arrived at the main 
features of the buckling problem, but the calculations are rather involved. 
Timoshenko assumed that the center line of the deflected beam as well as 
the initial shape is a half wave of a sine curve and arrived at a very 
simple solution. The restriction of the buckling mode to the symmetrical 
one, however, sometimes gives the critical buckling load manyfold too 
high in a certain range of arch rise. 

The buckling criterion used by the authors of references 1 to 5 
is the classical one which is based on the stability with respect to 

^After completion of the present work, it was learned that Hoff and 
Bruce (reference 6) treated a similar problem from the point of view of 
dynamic stability. Part of Hoff and Bruce’s work coincides almost 
identically with the present report. 
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infinitesimal displacements about the equilibrium positions. But 
Friedrichs, in reference 5, also calculated the buckling load on the 
basis of Tsien’s energy criterion, which is based on finite displace- 
ments. The energy criterion yields a buckling load much lower than 
that obtained from the classical criterion. It is not evident which 
of these two criteria corresponds to the real practical situation. 
Therefore in this paper, both criteria will be used and the results 
will be compared with experiments. 

This work was conducted at the California Institute of Technology 
under the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. 


SYMBOLS 


A cross-sectional area of beam 

a rise of arch 

E Young’s modulus 

F dead -weight load (in section ’’Buckling Load Based on Karman 

and Tsien’s Energy Criterion”) 

H axial compression at ends of beam 

H 0 initial thrust in beam 

I moment of inertia (or second moment) of cross section of beam 

K = tAeI^/aL^ 



L span of beam 

M bending moment; positive when it tends to put upper side of 

beam in compression 

M 0 bending moment due to lateral forces alone 
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Q 

q 

% 

t 

u 

V 

w 


shearing force in beam; positive when 



produces 


positive moment 

lateral pressure per unit length of beam; positive downward 
(in negative y-direction) 

characteristic lateral pressure per unit length of beam 
thickness of beam 
strain energy 

change in thrust in lateral support 

total load beam can sustain without buckling 


y, y 0 

a 

a' 

A 

5 

P„ 


actual and initial curve of center line of beam, respectively 
spring constant of arch support 
spring constant of lateral support 
distance spring-supported end of beam is displaced 
deviation ratio ^a m /a^ 


radius of circular -arc arch 
% bending stress in beam 

a c axial stress in beam 

a p = (j^Ei/l^/a 


0 total energy for dead-weight loading 

Subscripts: 

class classical criterion 

cone concentrated loading 

cr critical 
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energy energy criterion 

exp experimental 

max maximum 

sine sinusoidal loading 

unif uniform loading 

Nondimens ional coefficients: 
Let 


00 

r.-Z-* sin EH 


m=l 


-t 

m=l 


t^L sin 


mitx 


q = q^x) 


Then 


X 


m 


^4 

~ 



B 


m ~ 



R = 


qo^^ fX 

2jAei Vi 


HpL 2 

Jt 2 EI 


3 


a 


a + 


EA 

L 


= 2a 'L 3 
jAeI 


( m - l> 2 j 3> • • •) 
(m=l, 2, 3> • • •) 
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GENERAL ANALYSIS 


Consider a thin curved beam of small curvature, one end of which 
is hinged, that is, it is free to rotate but is fixed in position, while 
the other end of the beam is attached to a spring, with a spring con- 
stant a. When the spring-supported end is displaced by a distance A, 
the thrust induced in the spring will be 

H = H 0 + oA (1) 


where H 0 is the initial thrust built in the beam.^ (See fig. 3.) 
Before the application of the lateral load q(x), the axial load in the 
beam is H 0 and the beam center line is represented by the following 
Fourier series: 


y 0 


00 

> T • nutx 

l a-m-k Sln — f — 


m=l 


(2) 


as 


Under the lateral load q, the displaced center line can he written 


00 

y (3) 

m=l 

Assume that |y Q | and |y| are much smaller than L, and hence 
|a m | and [b^ are much smaller than lj that the beam is made of 
homogeneous material, of constant cross section, and with small curvature 

so that (dy/dx)^ is negligible in comparison with 1; and that the 
thickness of the beam is much smaller than the radius of curvature of the 
beam. Then the usual beam theory gives 



(4) 


No generality is lost by treating this case of one end spring 
instead of the case with both ends of the beam elastically supported 
because the springs at both ends can be replaced by a single spring at 
one end. 
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where M is the increase in bending moment due to the application of 
q(x). From statics, 

M = Q dx - (Hy - Hoy 0 ) (5) 

JO 

Substituting equation (5) into equation (4) and differentiating twice, 
remembering that ^ = -q and that the axial thrust in the beam H can 

be regarded as constant by the assumption of small curvature, the equa- 
tion of equilibrium is obtained: 


EI My-2 a> + H il 

dx 4 cbc 


Hr 


a 2 y 0 

dx 2 


= -q 


( 6 ) 


To find the thrust H, it is noted that the shortening of the 
center line of the beam is 



(7) 


where small quantities of higher orders are neglected. It is assumed 
that the end support spring is rather strong, so that A is very small 
compared with L. (Otherwise the problem becomes one of a simple bending, 
with no possible difficulty.) Hence 



( 8 ) 
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On the other hand, the deflection A is connected with the spring 
constant a by equation (l). Eliminating A between equations (l) 
and (8), substituting equations (2) and (3) for y and y Q into the 
result, and integrating, there is obtained 


H = H 0 + (3 


tt 2 EA 

4 


51 m2 ( a m 2 



( 9 ) 


where 


3 = 


a 


a + 


EA 


( 10 ) 


Substituting equations (2) and (3) again into equation (6) and using 
equation (9), there is obtained now the equation of equilibrium expressed 
in terms of the Fourier coefficients: 


-q = 



m (bm ' a m) 



- fcm) 


3 


jAeA 


5 ~ m2 (% 2 




2, . mnx 

m b m sin — — 

JL 


sin 


nnrx 

L 


(ID 


The boundary conditions are already satisfied. 
Expand q = q c f(x) into a Fourier series: 


where 


q = q^jfCx) 



, . mitx 

k m sln — 



sm 


mTtx 

L 


dx 


(12) 
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On equating the coefficients of the corresponding terms in the right- 
hand sides of equations (11) and (12), there is obtained a set of an 
infinite number of simultaneous equations: 


3 



n 





( b m ~ a m) - "V 


(m = 1, 2, 3, . . .) ( 13 ) 


To simplify the expressions, introduce the following notations: 




~2~ 





Then equations (13) become 



jt 2 EI 


(14) 


Bm 


00 00 
Y n2B n 2 - 3 YL “ 


X n 2 + m 2 


- S 


n=l 


n=l 


> 


- k m + ^m(m 2 - S) 
m^ 

(m = 1, 2, 3, • • .) (15) 


Here and ^ represent the rise of the arch, being half the ratio 

of the amplitude of the mth harmonic in the initial and the deflected 
curve to the radius of gyration of the beam cross section; R is a 
dimensionless quantity specifying the lateral loading; and S is the 
ratio of the initial axial compression to the Euler column buckling load 
of the beam. Now f(x), A. n , and S are known in the problem; it 
remains to find the relation between R and Bm, from which the corre- 
spondence between the load and deflection can be traced and the stability 
of the beam determined. 


Sometimes the Fourier series of the moment curve converges much 
faster than that of the loading itself. In such cases it is advantageous 
to use equations (4) and ( 5 ) directly instead of equation (6). Let the 
static bending moment of the lateral loading alone be written as M 0 : 
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M o = (16) 

where ^ is a characteristic lateral pressure with the physical dimen- 
sions of force per unit length of the heam. Let F(x) be expanded into 
a Fourier series, so that 


00 

M 0 = Qo l2 H K m sin SH 


m=l 


(17) 


where 


K m = f 


PL 


F(|) sin d£ 


Following the same reasoning as before, one arrives at the equations: 


Bmlp Y— n2 ®n 2 " P n 2 4 n 2 + 

\ n=l n=l 


m 


- S ) = -J^KjnR + X^m 2 - S) 


(m = 1, 2, 3, • . .) (18) 


B oth equations (15) and (18) will be used later. They are a system 
of an infinite number of simultaneous equations for which a general 
treatment is not known. However, there are many important cases where 
the number of equations can be reduced into a finite number; then a 
complete discussion is possible. Several examples will be given below. 

Equations (4) and (5) may be written as 

dfy H _ d2 y 0 M 0 (x) H q 
dx 2 EI 7 dx 2 + EI + El y ° 


= G(x) 


(19) 
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where G(x) is a known function. The general solution is 


y = C-l cos vx + C 2 sin vx + i I G(t) sin v(x - t) dt 

4 JO 


( 20 ) 


where 



The constants Cp and C 2 must he determined according to the boundary 
conditions at the ends y = 0 for x = 0 and L. The solution y(x) 
can then be substituted into equation (8) and 'v computed. This gives 
a relation between v and the external load. Biezeno and Friedrichs 
based their calculations on this relationship. Marguerre, on the other 
hand, used the energy principle and the methods of Ritz and Galerkin to 
obtain approximate solutions. The method of the present paper, based 
on the Fourier analysis, is due to the work of Y. S. Huang. 3 Recently, 
the same method was used by Hoff and Bruce (reference 6). 

It is clear from equation (20) that the deflection and the critical 
load are continuous functionals of y Q (x) and M 0 (x). Hence infini- 
tesimal changes in yo(x) and M 0 (x) would always cause an infini- 
tesimal change in the critical load. 


SINUSOIDAL ARCH UNDER SINUSOIDAL LOADING 


Consider the simplest case of a low sinusoidal arch subjected to 
a sinusoidal load distribution: 


y Q = a^L sin — 


q = qQ sin -± 


(0 < a 1 « 1) 

> 


( 21 ) 


J 


Professor of Aeronautics, Central University, Nanking, China. 
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The general equations of equilibrium ( 15 ) then become, in this 
particular case. 




3 i(p ]l n2 * 2 - ^’ 2 


n=l 


- 3Xq^ + 1 - S I = -R + Xq(l - S) 



£ 


n=l 


PXq 2 + 4 



( 22 ) 



n=l 


2n 2 


n B 


co 2 . J2. 

p A, ^ + m 



= 0 


This set of equations must be solved for B^ In order to get a 
qualitative investigation into the nature of the solution, first consider 
the simplest case of an arch rigidly hinged at both ends, so that a = <» 
or 3 = 1, and with zero initial axial thrust S = H Q = 0 . The more 
general case will be considered later. 

An obvious set of solutions of equations (22) is 
Bg = B^ = . . . = 0 

B l 3 - b i(£l 2 - 1) = \l - R 

If the relation between Bq and R is plotted, the curves in 
figure 4 are obtained. Depending upon the values of Xq, there are 
several possibilities: (1) If Xq ^ 1, the curves have monotonic slope; 

consequently, they determine the load-deflection curve uniquely. There 
is no question of instability. (2) If Xq > 1, then there are two real 
extremes and, for values of R in between these extremes, every loading 
may have three possible positions of equilibrium. Following, for example, 
curve I V in figure 4, the deformation of the beam can be traced as follows 
When the lateral loading is gradually increased from the starting point a 
the deflection gradually increases according to curve IV (the rise of the 


(23) 
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arch decreases). When the point M is reached, any further increase 
of loading will make the beam jump to the configuration corresponding 
to the point N and then follow the right-hand branch of the curve. 

In between M and N, any increase in deformation needs no addition 
of loading and therefore is unstable. Hence M is the critical point, 
with the critical condition given by 


dR 

dB! 


= 0 



< 0 


(24) 


From equations (23) and (24), the critical values of B]_ and R can 
be obtained: 


N 


Xi - l 


cr 




R cr = X 1 + 


#K x i 2 - x ) : 


(25) 


J 

If \j_ < 1 , R cr is imaginary; hence no instability will occur. 

This checks with the former discussion based on the uniqueness of the 
load-deflection curve. 


The above solution, equations (23), however, is not unique. Equa- 
tions (22) can allow a solution with one B n , in addition to Bj_, to be 

different from zero.* 1 In this case 


B l( B l 2 + n2s n 2 ) ‘ ( x i 2 - i)b x = X ± - R 
B ± 2 + n 2 B n 2 = A-! 2 - n 2 


(2 6 ) 


^These two cases exhaust the possibilities, as can be seen by writing 
down the rest of the set of equations (22), which gave the result that all 
other B's must vanish. 
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have the solution 





r (27) 


Equations (27) indicate that B n can exist (with real value) only in a 
definite range of R. The deformation history of a beam subjected to 
gradually increasing lateral loading can now be traced as in figure 5: 
Along ab, B n = 0, the curve is that of equations (23). Along be, 

B n ^ 0, the deflection curve becomes 


y = bpL sin -j- + b n L sin ^ (28) 

If the point b is real and lower than M, then it is the critical 
point where the beam will have a tendency to buckle. The point b is 
given by 



B n - 0 


y (29) 


R=kp-(n^-l) ^2 _ 


Equations (29) will yield the lowest critical value if the following 
conditions are satisfied: 

(1) R, Bp, and B n are real 

(2) The R given by equations (29) is less than the R given by 

equations (25) 

(3) The Bp given by equations (29) must be greater than that 

given by equations (25) ; otherwise, the beam will buckle in 
the first mode, at point M 
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(4) The particular number n is so chosen that the corresponding 
R cr is a minimum 

Conditions (l), (2), and (3) are satisfied if and only if 



(30) 


Condition (4) is satisfied only if n = 2. Hence the complete expression 
for the critical loading is obtained: 


The relation between the critical loading and the beam-rise ratio is 
illustrated in figure 6. The solid lines are the actual critical condi- 
tions. The dashed lines are either imaginary or not the lowest critical 
load. 

It is interesting to note here that for < {575 * 2.345 the 
buckling mode of a low sinusoidal arch is symmetrical but for X^ > 0-5 

the buckling mode imitates that of a high arch, for which the deforma- 
tion is essentially inextensional. As illustrated in figure 7, the arch 
deflects (flattens) at first under the increasing lateral loading from 
the initial position I to the state II, when the second mode Bg starts 
entering into the picture. The mode of the beam during buckling, when 
it jumps from the upper to the lower side, is a curve like III in 
figure 7- 


Still restricting this discussion to the simple case of a sinusoidal 
arch under a sinusoidal loading and with fixed hinged supports at both 
ends, let an initial compressive force H 0 act on the beam, so that 




(3D 


R C r - + 


3 /xp 2 - 4 


(Xi > 075) 


EFFECT OF INITIAL AXIAL COMPRESSION 
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is different from zero, S being the ratio of the initial axial compres- 
sion to the Euler column buckling load of the beam. The equation of 
equilibrium is given by equations (22) with 3=1. The solution of 
this set of equations is again either 


B-^ ^ 0, ~S>2 ~ — • • . — 0 


or 


Bp ^ 0, B n ^ 0, all the other B^s vanish 


One is led to the following conclusions: 


For Vl - S < X ± < \fW3 - S, 



and for X 1 > - S, 

R cr = (1 - S)Xp + 3/kp 2 + S - 4 


( 32 ) 


The effect of the initial axial compression is included in this formula. 
As expected, the increase of the initial axial compression will decrease 
the critical load, as can be easily verified by the fact that 


SR cr 

8S 


< 0 


( 33 ) 


for the full range of S, 0 ^ S ^ 1 (S cannot exceed 1). Furthermore, 
the lower limit for instability is now 

H = A - S (34) 


For X]_ smaller than this value, the bar is stable; no buckling is pos- 
sible. This lower limit decreases with increasing S until S = 1, when 
the beam will fail as a simple Euler column, R cr becoming zero. 
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The values of the critical load R cr as a function of X^, with 

values of S as parameters, are given in figure 6 and table I. A 
clearer presentation of the effect of S is a curve of the change in 
the critical load ^AR cr jg against Xp, where 

(^cr) s - ( R cr) g=0 - ( R cr) s , s (35) 

This is given in figure 8. 

From equations (32), it is seen that when Xp is large, say, with 
magnitude of the order of 2.5 or larger, ^AR cr ) s is almost linearly 
proportional to S. As a crude approximation, one may take 



INITIAL SHAPE OF ARCH OTHER THAN SINUSOIDAL 


In order to find the effect of the irregularities in the initial 
shape of the arch on the buckling load, some simple cases of low arches 
whose center lines are nonsinusoidal will be considered. By comparing 
such solutions with the previous one, the significance of such varia- 
tions in form can be estimated. Let the initial shape of the center 
line of the arch be given by the equation: 


y Q = apL sin — + a m L sin 


mffx 

L 


(37) 


(A few examples are shown in fig. 9*) Assume again for simplicity that 
the lateral loading q is sinusoidal, given by equation (20), and that 
the ends of the arch are hinged and without initial thrust, so that 
H q = S = 0 and p = 1. The fundamental equations (15) become 
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B l() n 2 B n 2 - - m 2 ^ 2 + lj = -R + A.-, 


Bml / n2fi n 2 - ^-l 2 - m 2 ^ 2 + m 2 ) = m 2 )^ 


( 38 ) 


bJ) ^n 2 * X l 2 - m2 \n 2 + k2 = 0 (for all k j 1, m) 


Again two possibilities exist: (l) A solution consists of f 0 and 

Bm i 0, but with all other B's vanishing; (2) a solution with one B k , 
other than B^ and B m , different from zero. They must be discussed 
separately. 

In the first case, Bm and R may be regarded as functions of B^ 
and the second of equations (38) differentiated to determine dBm/dBi. 
From the sign of dBm/dBi it can be observed that, when the load R is 
gradually increasing, the amplitude of Bm (i.e., | B m | j will increase 

irrespective of the initial sign of A^. Furthermore, by differentiating 
the first of equations (38) to obtain dR/cLBi, it can be observed that, 
in the prebuckling stage, the amplitude will decrease when the 

load R increases. Hence the critical condition is given by 



( 39 ) 


Carrying out the differentiation and reducing, the equation governing 
B m at the critical condition is obtained: 



+ cBni + d = 0 


(!» 0 ) 


where 


2X ± 2 + 2m 2 X m 2 - 3m 2 + 1 

2(m 2 - 1) 




31^ 

2(m 2 - 1 


c 
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Equation (4-0) can have at most two real roots. If the two real 
roots are different, then the one nearer to A, m Js the true critical 
value provided that the corresponding ( B l) cr and R cr are also real. 

If equation (40) has no real root, then there is no critical load and 
the beam is stable. 


With the critical value of B m so determined, the critical values 
of Bj and R can be obtained from equations (38) as follows: 


nj 


2 m V 


( B m) 


+ 


cr 


( B m) cr " V + 


3 


(4l) 


"cr ' M. * <” 2 - «fr)er - 


It is interesting to note here that the critical load is independent of 
the sign of A, m . This is so because a change in sign of A^ changes 
the sign of the roots (B m ) cr of equation (40). But since (B m ) cr /X m 
does not change sign, ^ from the first of equations (4l) is not 

affected by the change in sign of A^. Hence the conclusion follows 
from the second of equations (4l). This is rather unexpected. It shows 
that under sinusoidal loading the two apparently different curved beams 
in figures 9(b) and 9(c) have exactly the same critical load. 


Equation (40) can be solved graphically or numerically. The results 
of such calculations for the cases m = 2 and 3 are given in figures 10(a) 
and 10(b). The magnification of the amplitude of the higher harmonic, 
initially at A. m , into (B m ) cr at the critical point, is clearly seen 

from figure 10. The reduction of the critical load due to the presence 
of A. m will be discussed later when the second possible solution is 
obtained. The parameter used in the curves of figure 10 is not A^ 
but the deviation ratio: 


^ _ V _ 

x i a i 


(42) 


This ratio indicates the deviation from a sinusoidal form better than 
the parameter A^ itself. 
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It remains to discuss the second possible solution which includes 
one nonvanishing B k (k ^ 1, m). In this case the solution of equa- 
tions (38) is 


B 1 ‘ ,2 


R - Xi 


Bm = 


_ 1 


-m 2 X m 

v 2 2 

k - m 


k 2 B k 2 = X-l 2 + m 2 ^ 2 - k 2 - 


m°^n 2 (R - \ x ) 2 


(k 2 - m 2 j 2 " (k 2 - l) 2 


(43) 


The relation between B k and R is again an ellipse of a similar 
nature to that for a sinusoidal arch under sinusoidal loading. The 
instant when B k will appear is the critical point. Hence the condi- 
tion B k = 0 leads to 





iXj 2 - k 2 + m 2 ^ 2 


OJ 1 

% 

1 

CVJ 

X 

1 

i 1 


(44) 


This will lead to a fundamental critical value if the four conditions 
enumerated under equations (29) are satisfied. Whether equations (4l) 
or equation (44) gives the critical load depends on the initial shape 
of the beam. 

If m = 2, equation (44) always gives a higher R cr than equa- 
tions (4l) . Hence the critical load is determined by equations (4l). 

No B3, Bjp and so forth can appear during buckling. 

If m > 3, equation (44) with k = 2 gives the lowest R cr 
provided that is greater than a certain constant, say ( X^ ^ . For 

4]_ less than equations (4l) give the lowest R cr . The point 

( ^l)o is ^e point of tangency of the curves of R cr against X-^ 
computed according to equations (4l) and (44), respectively. 

Again it is evident from equation (44) that the critical load is 
independent of the sign of X m . 
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The combined results of equations (4l) and (44) are shown in 
figure 11, and the numerical results are given in tables II and III. 

In table II, (B m ) cr and R cr computed according to equations (40) 

and (4l) are listed. Comparing tables II and III, it is seen that in 
certain ranges of Xp, equations (4l) give the lower R cr , while in 
another range equation (44) gives the lower R cr . Furthermore, at 
smaller values of Xp, even if Xp > 1, (®l) cr and ^cr may become 

imaginary, as shown in table II. The physical meaning of this is that 
the process is then a continuous one. There is no sudden change of 
configurations. The beam, under bending, simply yields continuously to 
the increasing external load. 

These examples illustrate the serious nature of the effect of 
the Xjjj terms. When the loading is symmetrical, a very slight com- 
ponent of unsymmetry in the curved beam lowers the critical load con- 
siderably. For example, in case of a sinusoidal loading acting on a 
sinusoidal beam, an unsymmetrical second harmonic in the initial curve 
with an amplitude ratio of 1 percent in the initial form lowers the 
critical load by approximately 10 percent. The buckling mode is always 
unsymmetrical if the initial shape of the arch contains unsymmetrical 
modes. 

On the other hand, for a symmetric leading, the effect of higher 
harmonics that are symmetrical is much less pronounced. A similar 
effect should be expected when the beam is sinusoidal but the lateral 
loading deviates from a sinusoidal distribution. 

An important special arch form is a circular arc with a radius p Q . 
Within the present approximation, there may be written 



(45) 


This corresponds to an arch rise of 



at the center. The coef- 


ficients Aju form a rapidly decreasing sequence. In fact. 



a 2n = 0 



1 


(2n + l) 3 1 



^n+l " 
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The effect of the higher harmonics is negligible. If a^, a y, and so 

forth are neglected, then the R cr (sinusoidal loading) of a circular 
arch can be found from figure 11(b) or table III (m = 3) by taking 


The difference in R cr is readily seen to be small. 

To illustrate the fact that a^, ay, and so forth may be neglected 

without causing appreciable error, the case of the unsymmetrical buckling 
mode will be considered. Equations (43) should be modified, when k = 2 
into 




n 


n 


m(m^ - 4) 


Row 







m=3,5,7, ... 


m— 3 • • 


where 


*1 = 0.01436 


e 2 = 0.07325 
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The critical load is given by the equation: 

( 3 ~ ^ + t 1 + e 2 = i 1 + *l) x l 2 ~ k 

Neglecting the effect of X 5, X7, and so forth on R cr is to neglect 

the effect of and €2 on the root R of this equation. It is 

clearly justifiable. 


UNIFORMLY DISTRIBUTED PRESSURE 


In this section the critical load of a sinusoidal arch under 
uniformly distributed pressure will be discussed. From the results of 
the preceding sections, it is expected that the deflected curve of the 
arch would not remain sinusoidal and that an unsymmetrical component 
would in general enter into the buckling mode. For simplicity, again 
consider a simple sinusoidal arch, with ends hinged and without initial 
end thrust, so that 3 = 1 , Xg = X3 = . . . = 0 , and S = 0 . The 

lateral pressure is denoted by qQ per unit length of the span. Hence 
the bending moment in the beam, due to the lateral forces alone, with- 
out counting the contribution of the axial thrust, is 


M o = i <lo x ( L * x ) 


= < 3 o 


hlf 

n 3 


1 


n — i^ 3,5, 


1 . nflx 

n 3 Sin L 




It is convenient here to use equations (l8) because the Fourier series 
of M 0 converges much faster than that of the lateral loading itself. 
From equations (18), there can be obtained 


n 2 B n 2 - X x 2 + m 2 ) = - JL R + 8-^X-l 
\ ' JtnP 


1 y 3, 5 , • • •) 


m\ 


n^B n 2 - X-^ + m 2 j = 0 

(m = 2, 4 , 6, . . . ) 



Sim = 1 if m = lJ 

5 lm = 0 if * f 1 - 


( 47 ) 
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It is evident that when the load is applied, R ^ 0, all the B's 
with an odd subscript would in general differ from zero. It is also 
clear from the second of equations (47) that only one of the B’s with 


would be separately treated. 

Consider first the simpler case in which one of the B's with an 
even subscript is different from zero. In this case the deflection 
curve of the beam is unsymmetrical. Let the nonvanishing B be Bpk, 
where k is an integer. From the second of equations (47), then 


Obtained relating Bp^ with R. This relation is an ellipse, as in the 
section ’’Sinusoidal Arch under Sinusoidal Loading.” The critical condi- 
tion is reached when 


an even subscript can differ from zero, because 
equal to two different values of - m^. As 



cannot be 


As before, these two cases 



m 


Hence from the first of equations (U 7 ) , 




p 

Squaring Bm, multiplying by m , and summing,' there is obtained 




(50) 


O Q 

Equating this to - kVr according to equation (W), an equation is 


d^2k 
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which implies that 


B 2k = 0 


(51) 


With condition (51) , the critical load R cr is given by the following 
equation derived from equations (48) and (50): 


16 R 2 1 
** . L m*>(m 2 - I* s f 


8k L R 

n(l - 4k 2 ) 2 


1 



- 4k 2 
(52) 


p 

The series in the coefficient of R converges very fast. If all terms 
except the first one are neglected, the error is less than l/2 percent. 
Hence equation (52) is essentially equivalent to 


R - x l) 2 = - 1 ) 2 ( x i 2 - 4k2 ) (53) 


Comparing this equation with equations (29) for the case of a sinusoidal 
arch under sinusoidal loading, it can be seen that they are almost 

ij. 

identical except that R in equations (29) is replaced by — R and n 

is written here as 2k. One of the roots of equation (53) which would 
represent the critical load on the beam must satisfy the four conditions 
stated below equations (29). In a manner completely analogous to the 
treatment of sinusoidal loading under equations (29), one concludes that 
k must be equal to 1 and that the solution exists only when is 

equal to or greater than ^5- 5* The critical load is then 

R cr = T^l + 3^i 2 " ^ ) = ir( R cr) sine ( X 1 = l/^l) (5*0 


where (r.™^ . means the critical R of a sinusoidally distributed 
\ CI /sine 

lateral pressure. 

If the full series in the coefficient of R^ in equation (52) were 
taken, then, since k = 1, and 


d— 3,5,7, • . . 


■ 1 to = ^-977 x I0~ h 

m^(m 2 - 4) 
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equation (54) is modified by a factor of approximately (l - 0.005), 


or 


E =r - 0.995 £(R cr ) 


sme 


(*1 £ ^3) (55) 


Turning now to the other possible solution, that all the B's 
with an even subscript vanish, one sees by analogy to the case of a 
^dnusoidal arch under sinusoidal loading that this mode of deformation 
would lead to a critical buckling load only if \ 1 is sufficiently small. 
Let 


n=l ,3,5, . . 


2 P 

n% 2 = C 


(56) 


Then equations (47) give 


Bm = 


C - X.3 + m^ l jtm 


R + 5-^Xql 


(m = l, 3, 5, . . .) (57) 


From equations (57) compute m 2 ^ 2 and sum: 


■I 

m 


m 2 Bm 2 


_ 16 r 2 
' *2 R 


r 


8 

jr 


- R\ 1 


m=1 >3,5, ... 
1 


m 


- ^2 + m 2)2 

4l 2 


C ~ Xq 2 +1 C - Xq2 + q 


(58) 


is gives a relation between C and R but is rather useless because 
of its complexity. A more practical solution can be obtained by suc- 
cessive approximation. According to equations (57), for a given R 
Bm decreases rapidly with increasing m. As a first approximation’ 
then, neglect the effect of B3, B5, . . . and obtain from equa- 
tions (47), m = 1, the equation of equilibrium: 


B l 3 " ( X l 2 - l) B l = - — R + Xq 
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which is again almost identical with equations (22) for the case of a 
sinusoidal arch under sinusoidal loading, except that R in equa- 
tions (22) is now replaced by ^ R. Hence analogously. 



R 


cr 



( 59 ) 


= rMsi 


sme 


(is*-!* /5T3) 

J 


For the second approximation, neglect the effect of B^, By, and 
so forth, but consider B3. Now equations (47) may be written as 



Hence at the buckling point, 



- — - R, 
27 n 

_ 

" n R cr + X 1 + 


cr 


8 ( b i) 


cr 


x i + \j^i 2 ~ - 1 ) 3 

( l 8 H 2 - 23 l *)(Bl)cr 


Substituting into equations (60), which now become a relation 
combining B]_ and R, and using the criterion dR/dBp = 0 for buckling, 

there is obtained 


28 


MCA TN 2840 


( B l)cr * fe( k l 2 - 1 - “) 



^ (61) 


where 



—i2 

3 


(62) 


Since k is always positive, the critical load R cr given by equa- 
tions (6l) is always smaller than the first approximation given by 
equations (59) • But the difference is really very small because 
k « !• Values of k are given as a function of in table IV. 


vergence of the successive approximation is very good. From a com- 
parison of equations (62) and (59), there appears no need for further 
approximations . 

It can be concluded that, with an error less than l/2 percent, the 
critical value of R for a uniformly distributed loading is equal 
to ff/4 times that of a sinusoidal loading. 

Interpreting the result somewhat differently, compare the total 
load that an arch can carry when the load is distributed first uniformly 
and then sinusoidally. Let W be the total load. Then since 


Since (Bm) cr decreases very fast with increasing m, the con- 


W unif ~ 


and 



and since R cr is based on q^,, 



( 63 ) 
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Expressed in words, if W is the total lateral load an arch can 
sustain without buckling when the load is distributed uniformly over the 


load is distributed sinusoidally. Thus concentrating a load toward the 
center of the arch lowers the critical buckling load. 


The case of a concentrated load acting at the midpoint of the span 
can be analyzed in the same manner as for the case in the preceding 
section. Only a very brief explanation will be given below. 

Assume again that the arch is initially sinusoidal, rigidly hinged 
at both ends and without initial end thrust, so that (3 = 1, 

= . . . = 0, and S = 0. The lateral load is written as 


For the unsymmetric mode of buckling, if this mode is possible, the 
lowest R cr occurs when B 2 f 0, which implies that 


Q 

span, then the same arch can sustain only a total load of — — W if that 


CENTRAL CONCENTRATED LOAD ON A SINUSOIDAL ARCH 



(64) 


The bending moment due to the lateral forces alone is 





i Wx for 




m-1 




| W(L - x) for 


m— 1 , 3 > 5 j • • • 


(65) 


The equations of equilibrium are 


B, 


m 



- Xi 2 + 



(m - 1 , 2 , 3, ■ ■ •) (66) 



n 
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and the critical load is given by the smallest positive root of the 
equation: 


4r £ 


m=1 ,3,5, 


^ o - — R - | X.-, 2 + 4 = 0 

m 2( m 2 _ 4)2 9 9I 


Letting 


L 


1 + € 


m=l,3,5, , 


m< 


2( m 2 _ 4)2 9 


where € is approximately O.Ql+ 09 , 


or 



( 67 ) 


The numerical results of equation ( 67 ) which are used in the testing 
program are tabulated in table V and compared with ^( R cr) s i ne in fig- 
ure 12. For the symmetrical mode of buckling, steps analogous to those 
in the preceding section lead to 



(M. < ^53) (68) 


where 


k' = 9k 


k being the constant given by equation (62) and table IV. 
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Hence R cr for a concentrated load is approximately equal to one- 
half of R cr for a sinusoidally distributed load. 

As in the preceding section compare the load carrying capacity of 
a low arch with respect to various distributions of the loading. Thus 


(^cr)conc l/2 _ it 

( Wcr ) sine 2 / n ^ 

( Wcr ) sine . _8_ 
( Wcr )unif * 2 

( cr ) cone 8 x it _ 2 
( W cr)unif * 2 ^ " 


( 69 ) 


These ratios are within 2 percent of the corresponding ratios of 
the total loads causing equal center deflections of a simply supported 
beam under the three load distributions. This indicates that, for any 
symmetrical load distribution, the buckling load W cr is proportional 
to the total load (of the specified distribution) which causes unit 

center deflection of a straight simply supported beam.^ 


CENTRAL CONCENTRATED LOAD ON A N0NSINUS0IDAL ARCH 


Because the experiments to be described were carried out on a 
series of approximately sinusoidal arches with a central concentrated 
load, a more complete investigation of this case will be made. First 
the case f 0 and f 0 will be studied and the difference in 
the influence of A .3 on R cr for sinusoidal load and R cr for a 
central load will be shown. Then the case in which f 0, Xg ^ 0, 

and X 3 f 0 will be investigated. The results of the second case are 
more complicated and are used principally to show when the simple super- 
position of the effects of and ^3 is not possible. 


^This result was previously shown by Timoshenko (reference 1) for 
the case of symmetrical buckling mode. 
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For a pin-ended arch without initial thrust ( 0 = 1 , S = 0), the 
equation of equilibrium for a concentrated center load is 


Bm 


T~ n 2 (l n 2 - B n 2 ) - m 2 


2R . mir 2, 

= ^2 sin T “ m 


(m = 1 , 2 , 3 , 

If only Xq_ ^ 0 and 0, equations (70) become 


B l^l 2 + 9A-3 2 - n ^ B n 2 - = 


2R - X, 


Bgf A.^ + 9 \~ 


- n 2 B n 2 - k\ = 0 
n / 


B. 


s(^L 2 + 9 X 3 2 - Z! n 2 s n 2 - 9) = - | R - 9 X 3 


) (70) 


(71) 


For the case of buckling in the unsymmetric second mode one can 

V 2 2 

solve for ^ n B n from the second of equations (71). Substituting 

n" 

this into the other two of equations ( 71 ), solving the resulting equa- 
tions for B^ and B^, and again forming the sum n 2 B n 2 an equa- 

n 

tion is obtained connecting R, and B 2 . At the critical 

condition B 2 vanishes. Thus one arrives at an equation governing the 
critical load: 


4r £ 


I 


n— 1 , 3,5, . . 


n 2( n 2 _ k) 2 




sj ■ I X 1 2 * fy + ** * 0 
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Letting 


L 

n=l,3,5, • • • 


1 

n 2( n 2 - 4) 2 


k 

9 


k i 1.0409 


there is obtained 




X 3 


+ 


H - % m ) 2 * ^ ( 2 ^i 2 



(72) 


It is interesting to note that in contrast with the sinusoidal load 
case, the sign of A,^ I s important, and ^cr) conc no longer approaches 

^( R cr) S i ne unless ^ X ^ « Aq. In fact, the effect of X 3 on the 

ratio (Rcrj conc /(Rcrj sine can be appreciable. 

8 l 

The — A,^ terms in equation (72) arise from the cross product in 
the squaring of the right-hand side of the last equation of equations ( 71 ) 
to obtain B^. If* the case is considered in which only Aq and A^ 

are different from zero, there is no corresponding cross product and 
therefore it can be expected that A^ will have the same proportional 

effect on R cr for a centrally loaded arch as for a sinusoidally loaded 
arch. Physically this difference in the effect of and X .3 seems 

reasonable since the central load occurs at the maximum amplitude of 
the A .3 wave, but at a node of the A^ wave. 


For the case in which Aq, A, 2 , and A.^ differ from zero it is 

known from the section ” Initial Shape of Arch Other Than Sinusoidal” 
that buckling always occurs in the second mode and that the influence 
of the higher modes is small. Therefore all Bjq’s with m > 3 will 
be neglected. Letting p = 1 and S = 0, the equations of equilibrium 
are 


3 ^ 
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n 2 X n 2 - B x 2 - 4B 2 2 - 9B 3 2 - lj = 

B 2 ^L n2 ^n 2 - B i 2 " ^ B 2 2 - 9B 3 2 - 4 -J = -4X2 


2R - X, (73a) 


(73b) 


n2x n 2 ~ B l 2 - to 2 2 - 9B 3 2 - 9) = - | R - 9X, (73c) 


These equations are to be solved for the critical values of B^, 
b 2, B 3 , and R under the critical condition SR/SB-l = 0. The solution 

can be effected in the following steps: 

(l) Eliminate R between equations (73a) and (73c) and use 
equation (73b) to obtain an equation connecting B^ and B 2 : 



(2) Differentiate equations (73a) and (73b) with respect to B x 
and use the critical condition dR/cffi! = 0 and equation (73b) to obtain 

8B2/8B1 and 8B 3 /8Bi at the critical point. The results are expressed 
in terms of B^ B 2 , and B 3 « 

(3) These expressions for ^B 2 /dBi) cr and (dB^dBjJ are 

substituted into equation (73a) after differentiating it with respect 
to B^ • Using the critical condition dR/8Bp = 0 and eliminating B 3 

through equation (73b), an equation for (Bi) cr in terms of (B2) cr 

is obtained: 
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By plotting equations (74) and (75) a compatible solution can be found. 
This solution will not hold for Ag = 0, but it is valid for = 0 , 

although no simplification will result. The results for a series of 
arches with ^2/^1 = 0-005 and ^3/^1 = 0.04, which are representative 

of the test specimens to be described in the experimental section, are 
tabulated in table VI and plotted in figure 12 . Comparing this with 
figure 11 it is seen that the combined influence of Ag and A.3 is 
stronger than the sum of their separate influences for lower values of A.p 


can be used. This is not unexpected since for low values of X± the 

presence of Ag causes the mode of buckling to change from symmetrical 
to unsymmetrical and thus changes the influence of A^ on R cr . 

In figure 13 the process of loading is pictured for two example s in 
the above sequence of arches, one below the dividing value of Aq_ = ^5.5 
and one above. The changes in amplitudes of the three modes A.^ - b 1j 
B2 - Xq, and - X ^ are plotted as functions of the load R for 

= 2 and A,j_ = 4. It is to be noted that, for the lower value of 
A,]_, Bg does not increase rapidly until just before buckling occurs, 
while, for A,^ = 4, B 2 starts increasing rapidly at a point appreciably 
before the buckling point. For both cases B^ increases at an almost 
constant rate until just at the point of buckling. 


So far the ends of the arch have been considered as rigidly hinged. 
Since ideal rigid hinges cannot be realized in the testing machine, it 
is expected that some deviation in the experimental buckling load from 
the theoretical value may exist owing to the yielding of the supports. 

In order to obtain some quantitative measure of the effect of support 
displacement, an example of an arch with elastic supports will be 
considered. 

Assume that the supports are perfectly elastic. Let a be the 
spring constant of the support so that a displacement A would produce 
a thrust of magnitude ocA. Without loss of generality it can be assumed 
that one end is rigidly hinged, and the other is elastically supported, 
as shown in figure 3. The effect of the support rigidity on the 


but for higher values of 



the principle of superposition 


ELASTIC SUPPORTS AT ENDS 
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equilibrium is expressed by the parameter p, defined by equation ( 10 ). 
The equations of equilibrium are either equations ( 15 ) or (18). 

As an example , consider a sinusoidal arch loaded laterally by a 
sinusoidally distributed load of intensity q per unit length: 


y Q = a^L sin — 
q = qo sin ig 


The equation of equilibrium is given by equations (22). The solution 
obtainable in the same manner as in the section "Sinusoidal Arch under 
Sinusoidal Loading" is 


E cr = (1 - S)X X * ~ 1 * ^ 

(w <hi (Xl) °) 

. s ,/pX-, 2 - 4 + S 

^cr ~ ~ S)X.j_ + 3u p 

( X 1 > Mo) 


where (\ 1 J 0 is the smallest positive real root of the equation: 


( 76 ) 


(<% 2 - 4 + S ) = 2 F 3 '( PX 1 2 - 1 + S ) 3 


( 77 ) 


The effect of the nonrigidity of the support (p < 1 ) is shown in fig- 
ure 14 and table VII. The values of as a function of P are 

also given in that figure and table. The limiting case, a > » and 

P checks with the results of the sections "Sinusoidal Arch under 

Sinusoidal Loading" and "Effect of Initial Axial Compression." 

If the support offers no resistance to the axial thrust, that is, 
it is perfectly flexible, then a = 0 and p = 0. In this case there 
is no critical buckling load} the arch deflects continuously because the 

lower limit of A-i, 1 //P, below which no buckling can occur, now tends 
to infinity. 
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Similarly other loading conditions may be treated. For example, 
if 3 differs from 1, the ratios of R cr for a uniformly distributed 
load, a sinusoidally distributed load, and a concentrated load at the 
center are again, respectively, ir/4, 1, and l/2. 


LATERAL ELASTIC SUPPORTS 


In application to certain wing design problems, it is desired to 
investigate the effect of lateral elastic supports on the buckling load 
of the arch. As an example, consider an arch having an elastic support 
at the center, as shown in figure 3. Let a' be the spring constant 
of the support. Then the change in thrust in the spring is 

V = o' ( a* - (A’) 0 ] (78) 

where A' - (A') 0 is the change in the deflection at the midspan. No 
generality is lost by assuming (A') q to be zero, if initial thrust in 
the spring is counted as a lateral force. 

Now when the deflection curve of the arch is given by equations (2) 
and (3), 


L m-1 

(~l) ^ ( a m - fom)L 

m=l, 3 y 5 > • • • 


(79) 


The moment contributed by V is then (cf. equation (65)) 


m-1 


(m \ 2VL 

\ M o/V = - — 


ni— 1 , 3 y 5 y • • 


(.1) * _L Sin 


( 80 ) 


Combining equations ( 78 ), (79 ) > and (80) and adding (M q ) v to the right- 

hand side of equation (4), there is obtained, after some reduction, the 
equations of equilibrium (equations (l8) modified): 
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^0 n2s n 2 “ n 2 X n 2 + m 2 - sj = -ir^R + xjm 2 - s) + 


sm 


— JL 

2 2 
m 


n =l,3,5, • 


n-1 

2 


(-1) ^X, n - B n ) 


(81) 


where K m is given by equation (17) and 


-■Si 


Since a simply supported beam with a unit concentrated load at its 

3 3 

center has a deflection of -L L. ~ 2 L under the load ^ ig 

4o "J- h ii* I 

jt 

approximately the ratio of a' to the spring constant of a simply sup- 
ported beam having the linear dimensions of the arch. 

Consider a sinusoidal arch subject to sinusoidal loading. For 
simplicity let the initial thrust be zero and let the end hinges be 
rigid. Then if m ^ 1, 


3 = 1 

S = 0 

V = o 



K m = 0 
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The governing equations are 


n2s n 2 * + “ 2 J = + 5^! + 


r 


m-1 


(-1) * 4 n 


\-\ - 


L 


n=l,3,.. . 


n-1 

(-1)" 2 ” B n 


(if m is odd) 


(if m is even) 


(83) 


Let 


Q = - 

n=l,3, . . . 

p = Yi n2 v 

n=l,3,... 


n-1 

2 


(-D ^ B n 


(84) 


then 


Bi = 


- R 4- 2M-Q 
P + 1 - X.-, 2 


m-1 


B m - (-1) 


2HQ 


+ m 2 - X.^) 


(m odd and > 3) 


m* 


With these values for there is obtained 


Q = + 


R - X... 


2pQ 


P + 1 - \ 


± 4- m2 ( p2 + m2 - x i 2 ) 


( 85 ) 
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First calculate the critical load for unsymmetric buckling where 
B 2n t 0 for some n - Then according to equations (83), 


P = X-, 2 - 4n 2 


( 86 ) 


But P is also given by 

p - H An 2 

m 


. wv y 

m=l >3 , • • 


i 2 (m 2 - 4n 2 )‘ 


pl 


m 


R)(X! - B + tpQ) 2 

(l - kr?) 2 2n 


(87) 


Neglect all except the m = 1 term of the series and substitute 
the value of P from equation (86): 


- R + 2HQ.) 2 

(l - 4n 2 ) 2 


- 4n 2 - ^n 2 B 2n 2 


(88) 


The critical condition is again B 2n = 0. Solving equation (85) for Q 
and substituting together with B2 n = 0 into equation (88), one obtains 

fs : i,„2 + su h - ^ 

(4n 2 - l)(l + 2|ia) + 2 h 

where 


R C r - ^2. + “ l) 


a 



1 

m^(m^ - kn^) 
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This expression is a minimum for n = 1, hence 


a = -0.3087 

R cr ( Rcr ) sine + 3 + O.l^Tn^ 1 ” ( Rcr ) sine] 




(89) 


These values are given in table VIII and plotted in figure 15 for 
P = 0, 0.5> 1, 1 .5, 2, and 3* 


In this solution the effect of the higher modes (m ^ 3) on the 

force exerted by the spring, which enters by the series > — - , 

m 2 (m 2 - k) 

is included, but the effect of the higher modes in lowering the buckling 
load, which enters by the faster converging series \ i , is 

Z _ m 2( m 2 _ 4) 2 

neglected. In the analogous case of the arch with a concentrated central 
load this result's in a maximum error of 3 percent for < 10 and for 
this case it should be no more. 


Consider next the case of symmetrical buckling which occurs for the 
smaller values of X.^. As a first approximation neglect the effects of 
all the B^s except B^. Then from equations (83), there is obtained 
under the critical condition dR/dBp = 0 the critical load: 


R cr = Xj_(i + 2p) 


+ 





(90) 


A procedure similar to that used in the section "Initial Shape of 
Arch Other Than Sinusoidal" can be applied to find further approxima- 
tions. The results of such a calculation, with the effects of Bj_ and 

B3 included, are given in table VIII and are plotted in figure 15 . 
BUCKLING LOAD BASED ON KARMAN AND TSIEN'S ENERGY CRITERION 


It is well-known that the classical buckling criterion, on which 
the calculations of the preceding sections are based, leads to erroneous 
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results for cylindrical and spherical shells; while a fundamentally 
different criterion, first proposed by Karman and Tsien, whose latest 
version is given in reference J, gives much closer agreement with exper- 
iments. The criterion of K^rm^n and Tsien (henceforth referred to as 
energy criterion) is that the buckling load is reached when the total 
energy in a possible (buckled) equilibrium state is equal to the total 
energy in the unbuckled state. In other words, if the total potential 
energy is such that it is permissible for the structure to jump from 
the unbuckled state to a buckled state, then the structure will actually 
jump. 

Both the classical and the energy criteria have been applied to 
curved beams and shells. In some cases the classical criterion gives 
closer agreement with experiments; in others, the energy criterion gives 
better results. The reason, as pointed out by Tsien, is that in some 
cases the energy "hump" between two equilibrium states (one buckled and 
one unbuckled) of the same energy level is large and in other cases it 
is small. If the hump is small, the ever present small disturbances 
will enable the structure to jump from the unbuckled state to the more 
stable buckled state. Otherwise, this jump will be hindered. The 
crucial decision of the proper criterion depends much on what one means 
by a practical" experimental setup or a "practical" service condition 
of the structure . 

The energy criterion has been applied to the low arch problem by 
Friedrichs (reference 5) who found a great reduction in R cr based on 
the energy criterion from that based on the classical criterion. In 
order to decide which criterion actually applies to the buckling of low 
arches, the experimental setup to be described in the next section 
will be accepted as practical and the theoretical results will be com- 
pared with experiments. 

In applying the energy criterion, one must distinguish a constant 
deflection loading (a rigid testing machine) from a constant force (dead- 
weight) load. In the former case the change in total energy in buckling 
is just equal to the change in the internal strain energy, while in the 
latter case it is equal to the change in the strain energy minus the 
force times the displacement. However, a laterally loaded arch cannot 
buckle if the point of loading is not allowed to jump; hence only the 
dead-weight loading case will be considered. 

For dead-weight loading the total energy is 


0 = u - w (91) 

where U is the strain energy and W is the work done by the lateral 
loading. The energy 0 can be expressed as a function of the 
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deflection 5. Then according to the energy criterion, buckling would 
occur under a dead weight F provided that 


0( 8 l) - 0( 8 2) 
F ( 8 l) * F ( 8 2) 


(92) 


where 5^_ and 5 2 are two deflection configurations. Now the strain 

energy U, under the assumptions of the section "General Analysis," is 
given by 


PL 


U = - El 
2 


d 2 y 


idx dx 


ii °\ 2 dx + (" 2 - h ° 2 ) l 

2 / 2AE 


(93) 


From equations (2), ( 3 ), and (9), equation (93) becomes 


U = K m h j(Xm - B m ) 2 [l + §( A-m + B^ 2 ] + S (\ m 2 - B m 2 )l (94) 
m v L_ _J J 


where 


K = 



The work done by the external load in the buckling process is 


W 



q(y Q - y) dx 


For a sinusoidally distributed loading, 


( 95 ) 


W = 2KR(X. 1 - B ± ) 


(96) 
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while for a concentrated load at the midspan 


W = 4kr 


I 

m=l,3,5, . . . 


m-1 

(-1) 2 (A. m - Em) 


(97) 


The buckling load according to the energy criterion can then be obtained 
easily. 


Sinusoidally Distributed Loading on a Sinusoidal Arch 

It was shown in the section "Sinusoidal Arch under Sinusoidal 
Loading" that the only equilibrium position of a sinusoidal arch under 
sinusoidally distributed load is the one for which all the Bm's 
(m = 2, 3> • • • ) vanish. Hence if S = 0 (zero initial thrust). 



(X-l - Bp) 2 1 + + BJ 2 ] - 2R(\p - Bp) 


1 


> 

R = X-p + Bp(x.p 2 - l) - Bp3 


(98) 


The buckling conditions that 0^Bp') = 0(Bp") and R^Bp') = R^Bp") are 
fulfilled when R = X.^ at which 


| = Bl 2 - ^2 + |(Xp 2 - Bp2) 2 ( 99 ) 

and = 0^-B-J. A substitution of R = X 1 into the second of 

equations ( 98 ) gives the arch rise at the critical condition: 


or 


Hence 




- 1 


( B l) cr - 0 
^cr = ^1 


( 100 ) 
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Central Concentrated Load on a Sinusoidal Arch 
Assuming no initial thrust (S =0), from equations (94) and (97) 


K ~ ( X 1 ' B l) 2 £ + |( X 1 + B l) 2 




From equations (66) 



m-1 

+ 4R(-1) 2 


2R = X ± 




( 101 ) 


( 102 ) 


If all the Bja’s except are neglected, the above equations 

become identical with those for the sinusoidal loading if R is replaced 
by 2 R. Thus approximately, R cr for the concentrated center load is 
one-half of that for the sinusoidal load. This is the same approximate 
ratio as for R cr of the sinusoidal and the concentrated loadings based 
on the classical criterion. 


The ratio of R cr based on the energy criterion to that based on 
the classical criterion is plotted in figure l6 for sinusoidal loadings 
on a sinusoidal arch. This same ratio holds approximately for the 
central load on sinusoidal arches. 


EXPERIMENTS 


A series of pin-ended arches having rigid simple supports were 
loaded with a central concentrated load in the testing apparatus shown 
in figures 17 and 18. The ideal end conditions were approximated as 
closely as possible by supporting the arches on knife edges mounted in 
a heavy steel frame having a stiffness approximately 100 times that of 
the specimen. Allowing a 20-percent reduction in this stiffness due to 
the flexibility of the knife edges and fittings results in a value of 3 
equal to 0.988. A reference to the section "Elastic Supports at Ends” 
and figure lk shows that a maximum error of about 1 percent will result 
from considering the supports as perfectly rigid. 
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The knife-edge fittings were provided with sockets which alined 
the ends of the specimen with the knife edges. (See fig. 19*) 

The most critical problem in setting up the specimens for testing 
was the spacing of the supports. A looseness or an initial compression 
results in a change in the initial arch shape and an appreciable error 
in the buckling load. In the tests the spacing adjustment was made by 
a wedge controlled by a screw which was rotated until the play between 
the specimen and the knife edges was just eliminated. 

The specimens were cut from 24S-T3 and 75 S-T 6 sheets and milled to 
l/2-inch width. The strips were then rolled to the desired curvature 
on a three-roll roller. To reduce the effect of roll eccentricity 
several passes were made at each setting of the rolls, the rolls being 
indexed to a new position at the start of each pass. 

The curvature of each specimen was measured at 12 stations by a 
dial gage which could be read up to ten-thousandths of an inch, placed 
between knife edges 2 inches apart. These curvatures were numerically 
integrated to find the shape of the specimen for which a 12-term Fourier 
expansion (half -range sine series) was made. The first three coefficients 
are given in table IX. As a check on the accuracy of the method the 
central rise of the arch as predicted by the numerical integration was 
compared with the actual rise as measured with a vernier height gage. 

The difference was no more than 4 percent of the arch rise for each 
specimen measured. The central arch rise as predicted by the Fourier 
coefficients agreed with the numerical integration within 1 percent. 

The Fourier coefficients Aq, A^, and A, 3 were used in calculating 

the theoretical critical load. In such calculations use is made of the 
fact noted in the section "Central Concentrated Load on a Nonsinusoidal 
Arch" that, whereas for smaller Aq ( sa y> Xq < 2.4) the joint effect 
of X -2 and A.^ on R cr is not equal to the sum of the effects of A^ 
and A>3 separately, for larger Aq (say, Aq > 2.4) the effects of Xq 
and A,^ are superposable.. Hence for < 2.4 the more exact method 
of the aforementioned section was used, but for Aq > 2.4 the effects 
of A <2 and A^ were calculated separately and added together algebra- 
ically. The effect of A .3 is given by equation (72). That of Xg, 

according to the previous argument, can be obtained, percentagewise, 
from figure 11(a) or table III. 

Although no attempt was made to determine the arch shape during 
the loading process, visual observation showed that the test performance 
at least approximately agreed with the theoretical predictions. The 
gradual increase in the third mode with the load, resulting in a flat- 
tening of the arch and then a reversal of curvature for the higher values 
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of Aq, was noted. For values of Aq > 2.4 the rapid increase in the 

unsymmetric second mode just before buckling was quite evident. The 
clearest indication of the onset of buckling, however, was obtained by 
noticing the vibration of the specimen as the individual weights were 
applied. Even a very careful application resulted in a slight vibration 
in the fundamental mode. When the load approached within a few pounds 
of the critical load there was a rapid decrease in the frequency of this 
vibration. Further load applications were made in extremely small 
increments. 


The theoretical and experimental results are listed in table IX 
and plotted in figure 1 6 . In figure 1 6 , the ordinate is the ratio of 
R C r determined by the test to that computed theoretically according to 
the classical criterion. In the same figure, the dashed line shows the 
ratio of R cr given by the energy criterion to that given by the clas- 
sical criterion. This curve is based on the simple sinusoidal arch 
(^2 = X3 = o). For arches used in the experiment \r> and X ^ were so 


small that the variation of the ratio /(R™) does not 

\ / energy/ V /class 

vary much from the dashed curve of the figure. 


It is seen that the test results agree quite well with results 
based on the classical criterion for higher values of Aq but drop 
appreciably below them for the lower values. All the test values, 
however,, lie above the energy criterion curve. Although calculations 
for the series of arches representing the test specimens indicate that 
buckling would occur for X ^ 1.05, no buckling was observed for 

X 1 ^ 1.38. 


A calculation of the stresses in the specimens at buckling was 
made to determine if yielding occurred. With H 0 = 0 and 3=1 the 
axial compressive force is given by equation (9). Using the nondimen- 
sional notation it becomes 



■ TTp- I! m2 (v 2 - V 1 ) (103) 

AL, m 


For a sinusoidal arch with a sinusoidal load all the Bm's except 
B-l are zero at the critical buckling load and 2 = j(\l 2 " ^ 

for 1 ^ X x 2 ^ 5 .5 and (bJ ^ 2 = X-) 2 - 4 for X x 2 £ 5 . 5 . Therefore 
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| (l * 2^ 2 ) 

4 


(for 1 4 x 2 < 5 . 5 ) 


(104) 

(for A.-J 2 ^ 5 - 5 ) 


where P = jr^El/L^ is the Euler buckling load of the beam and ap = P/A. 
Thus it can be seen that the critical compressive force is just equal to 
the Euler load if = 1. As increases the critical force increases 

until it reaches the Euler load for buckling in the second mode. At 
this point the arch buckles unsymmetrically and the critical compressive 
stress remains constant for all higher values of A,j_. This performance 
is also typical of symmetrical arches with a central concentrated load, 
but for arches with a slight asymmetry, as is the case for the specimens 
tested, the value H cr /P = 4 is approached only as becomes large. 

The values of H cr /P for a series of arches are given in table VI. 

The maximum bending stress at any point x is given by 



(105) 


where t is the thickness of the specimen. In terms of the nondimen- 
sional Fourier coefficients this becomes 


a b 


rt^Et 

l 2 



“ 2 (*m - Bm) 


sin 


mTtx 

~ 


(106) 


The bending stresses at the midspan were calculated for the series of 
arches with A .2 = 0.005Xq. an( ^- ^3 = 0. 0*4-0 A.-^ which are representative 

of the actual test specimens. The results are shown in table VI together 
with the total maximum stress for t = 0.25. The total stress for any 
other thickness is obtained by multiplying the last column of table VI 
by the factor l 6 t^. 

All the specimens tested had maximum stresses well below the yield 
stress of the material at the buckling point. Yielding occurred in the 
post-buckling stage for all the specimens except those having the very 
lowest values of A.]_. 
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CONCLUSIONS 


A Fourier analysis has been used to solve the problem of buckling 
of low arches under a lateral loading acting toward the center of 
curvature. The conclusions may be summarized as follows: 

1. For a sinusoidal arch under a sinusoidal loading, the analysis 
gives a very simple exact solution for the nonlinear equation of equi- 
librium. The critical load can be expressed as a simple function of 
the beam dimension parameters. On the basis of the classical buckling 
criterion, it is shown that the buckling mode is symmetrical for arches 
having a nondimens ional parameter Aq less than ^5-3 and is unsym- 
metrical for Aq greater than V3- 5* This dividing value is affected 
somewhat by the initial thrust in the arch and the elasticity of the 
support. 


2. For arch shapes other than sinusoidal but under sinusoidal 
loading, it is shown that symmetrical deviations have only minor effects 
on the buckling load, while unsymmetrical modes of deviation cause 
serious reduction of the buckling load. The buckling mode is always 
unsymmetrical if the initial shape of the arch contains unsymmetrical 
modes. For sinusoidal loading the critical load is independent of the 
sign of A. m (m > l); thus a pair of different arches can have the same 
critical load. 

3. For a load distribution that deviates from sinusoidal, the 
unsymmetrical components again have serious effects. The critical load 
will be dependent upon the sign of A^m > 1). For symmetrical load 
distributions, the buckling loads are approximately proportional to the 
total loads (under the respective distributions) that are required to 
produce a unit deflection at the center of a straight simply supported 
beam without axial restraint. 

4. Comparison with experiments shows that the classical criterion 
of buckling is applicable for larger values of Aq, say, Aq > 3* But 
the classical criterion overestimates the buckling load for very flat 
arches. The experimental buckling load is always higher than that 
estimated according to the energy criterion of Karman and Tsien but has 
a tendency to approach that criterion as Aq decreases. For Aq > 1 

(with exact value depending on the initial thrust and support conditions), 
the arch deflects continuously and there is no buckling phenomenon.. 

California Institute of Technology 

Pasadena, Calif., January 24, 1952 
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TABLE I 


VALUES OF R cr AS A FUNCTION OF INITIAL THRUST AND ARCH RISE 


s \ 

(1 - S) 1 / 2 

1.0 

1.2 

1.4 

1.6 

1.8 

2.0 

0 

1.000000 

1.000000 

1.312338 

1.762039 

2.349955 

3.090387 

4.000000 

.1 

.853815 

. 912172 

1.232735 

1.680056 

2.263209 

2.997755 

3.900829 

.25 

.649519 

.798113 

1.120608 

1.562302 

2.137273 

2.862332 

3.755138 

• 50 

.353553 

.636083 

.950784 

1.379012 

1.938017 

2.645717 

3.520288 

1.00 

0 

. 384900 

.665108 

1.056166 

1.576551 

2.244738 

3.079201 

\*1 
s \ 

2.2 

(5.5 - s) 1 / 2 

2.4 

2.6 

3 

3-5 

4 

0 

5.096309 

6.019436 

6.379947 

7.583975 

9.708204 

12.116843 

14.392306 

.1 

4.990180 

5.765646 

6.251454 

7.413459 

9.474954 

II.818911 

14.035515 

.25 

4.833707 

5.392701 

6.053234 

7.154805 

9.123864 

11.371428 

13.500000 

.50 

4.580028 

4.792269 

5.709987 

6.716641 

8.535621 

10.624117 

12.606599 

1.00 

4.O69398 

3.674235 

4.983975 

5.817216 

7-348470 

9.124145 

IO.816653 
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TABLE II 


VALUES OF (Bm) cr , 


(Bi) cr , AND R cr COMPUTED FROM EQUATIONS (1(0) AND (1(1) 


(a) m = 2 


X 1 

*2 

Her 

(®2_)cr/^2 

( B l)cr 

^cr 

*2 

Her 

Her/ ^2 

Her 

R cr 


X.2/X1 = 0.01 

Xg/Xi = 0.05 

1.0 

0.010 

0.01334 

1.334 

Imag.l 

Imag. 

0.05 

O.06655 

1.331 

Imag. 

Imag. 

1.2 

.012 

.01773 

1.478 

0.3828 

1.312 

.06 

.08824 

1.471 

0.3785 

1.306 

1.4 

.011( 

.02372 

1.694 

.5656 

1.761 

.07 

.1170 

1.672 

.5635 

1.742 

1.6 

.016 

.03259 

2-037 

.7214 

2.348 

.08 

.1569 

1.961 

.7257 

2.297 

1.8 

.018 

.0471(4 

2.636 

.8660 

3.084 

.09 

.2121 

2.357 

.8886 

2.958 

2.0 

.020 

.07663 

3.832 

1.0110 

3-978 

.10 

.2817 

2.817 

1.0690 

3.689 

2.2 

.022 

• 1339 

6.088 

1.1947 

4.999 

.11 

• 3568 

3.243 

1.2699 

4.443 

2.4 

.02k 

.1990 

8.294 

1.4443 

6.036 

.12 

.4302 

3.585 

1.4808 

5.190 

2.6 

.026 

.2534 

9.748 

1.7077 

7.022 

.13 

.5002 

3-847 

1.6931 

5.919 

3.0 

.030 

• 3428 

11.43 

2.2099 

8.856 

.15 

.6307 

4.204 

2.110 

7.322 

3-5 

.035 

.4386 

12.53 

2.7936 

10.99 

.175 

.7840 

4.480 

2.609 

8.997 

4.0 

,040 

.5267 

13.17 

3.3467 

13.02 

.20 

• 9^71 

4.636 

3.096 

10.616 


X2/X1 = 0.1 

Xg/Xi = 0.2 

1.0 

0.10 

O.1325 

1.325 

Imag- 

Imag. 

0.20 

0.2606 

1.303 

Imag. 

Imag. 

1.2 

.12 

.1742 

1.451 

0.3640 

1.289 

.24 

.3349 

1.396 

0.2964 

1.240 

1.4 

.14 

.2264 

1.617 

• 5540 

1.692 

.28 

.4187 

1.495 

.4976 

1.562 

1.6 

.16 

.2910 

1.818 

.7235 

2.179 

.32 

.5094 

1.592 

.6668 

1.925 

1.8 

.18 

.3662 

2.035 

.8940 

2.724 

.36 

.6039 

1.677 

.8272 

2.309 

2.0 

.20 

.4473 

2.237 

1.0715 

3.298 

.ko 

.6998 

1.749 

.9838 

2.702 

2.2 

.22 

.5296 

2.407 

1.2542 

3.879 

.44 

.7955 

1.808 

1.1381 

3.096 

2.4 

.24 

.6108 

2.545 

1.4388 

4.455 

.48 

.8906 

1.855 

1.2902 

3.489 

2.6 

.26 

.6900 

2.654 

1.6226 

5.022 

.52 

.9848 

1.894 

1.4403 

3-879 

3.0 

.30 

.8432 

2.811 

1.9848 

6.130 

.60 

1.170 

1.951 

1.7354 

4.648 

3-5 

.35 

I.027 

2.934 

2.426 

7.470 

.70 

1.398 

1.997 

2.096 

5.590 

4.0 

.40 

1.205 

3.013 

2.816 

8.709 

.80 

1.622 

2.028 

2.450 

6.518 



X^/Xi 

= 0.3 



x^y = 0.4 

1.0 

0.30 

0.3829 

1.276 

Imag. 

Imag. 

0.40 

0.5002 

1.250 

Imag. 

Imag. 

1.2 

.36 

.4824 

1.340 

0.1111 

1.202 

.48 

.6212 

1.294 

Imag. 

Imag. 

1.4 

.42 

.5879 

1.400 

• 3746 

1.453 

• 56 

.7456 

1.332 

Imag. 

Imag. 

1.6 

.48 

.6968 

1.452 

.5429 

1.733 

.64 

.8716 

1.362 

0.3177 

1.612 

1.8 

.5^ 

.8071 

1.495 

.6907 

2.024 

• 72 

.9980 

1.386 

.4637 

1.853 

2.0 

.60 

.9208 

1-535 

.8093 

2.319 

.80 

1.124 

1.405 

.5916 

2.091 

2.2 

.66 

1.028 

1-557 

.9620 

2.615 

.88 

1.250 

1.421 

.7082 

2.331 

2.4 

.72 

1.137 

1.580 

1.0911 

2.911 

.96 

1.376 

1.433 

.8181 

2.571 

2.6 

.78 

1.246 

1.598 

1.2174 

3-205 

1.04 

1.500 

1.443 

.9236 

2.810 

3.0 

.90 

1.463 

1.625 

1.4640 

3.789 

1.20 

1.749 

1.458 

1.1242 

3.288 

3.5 

1.05 

1.730 

1.648 

1.7640 

4.510 

1.40 

2.058 

1.470 

1.3688 

3.881 

4.0 

1.20 

1.996 

1.663 

2.0581 

5.224 

1.60 

2.365 

1.478 

1.6048 

4.471 


•*-Imaginary . 
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TABLE II 

VALUES OF (Bmjcj., (Bi) cr , AMD R cr COMPUTED FROM EQUATIONS (1+0) AND (1+1) - Concluded 

(b) m = 3 


( B 3) c 


( B 3)cr/> 


( B l) 


cr 


( b 3) 


cr 


( B 3j 


cr/ x 3 


( B l)c 


X3/X1 = 0.01 


1.0 

0.010 

0.01125 

1.125 

Imag. 

1.2 

.012 

.011(02 

1.168 

0.3779 

1.4 

.014 

.01712 

1.223 

.561+7 

1.6 

.016 

.02069 

1.293 

• 7342 

1.8 

.018 

.021+90 

1.383 

.8622 

2.0 

.020 

.02998 

1.1+99 

1.0000 

2.2 

.022 

.0361+2 

1.655 

1.1265 

2.4 

3.0 

.024 

.01+1+60 

1.858 

1.2611 

4.0 

— 




— 


X3/X1 = 


0.05 


Imag. 

1.312 

1.762 

2.361+ 

3.088 

3.996 

5.083 

6.381 


0.05 

0.05623 

1.125 

Imag. 

Imag. 

.06 

.07001 

1.167 

0.3762 

1.308 

.07 

.0851+0 

1.220 

.5617 

1-750 

.08 

.1030 

1.287 

.7183 

2.323 

.09 

.1236 

1.374 

.8530 

3.035 

.10 

.1473 

1.473 

1.0025 

3.894 

.11 

.1758 

1.598 

1.11+25 

4.903 

.12 

.2098 

1.748 

1.2817 

6.054 

— 









X3/A1 = 0.1 


x 3/ x l 


A.+ = 0.2 


1.0 

0.10 

0.1112 

1.122 

0.11+63 

0.9969 

0.20 

0.2234 

1.117 

Imag. 

Imag. 

1.2 

.12 

.1395 

1.162 

.3693 

1.295 

.24 

.2760 

1.150 

0.3172 

1.255 

1.4 

.14 

.1696 

1.211 

• 5538 

1.716 

.28 

• 3321 

1.186 

.5110 

1.610 

1.6 

.16 

.2032 

1.270 

.7109 

2.21+9 

.32 

.3919 

1.225 

.6694 

2.036 

1.8 

.18 

.2408 

1-338 

.8587 

2.893 

.36 

.1+51+9 

1.261+ 

.8163 

2.516 

2.0 

.20 

.2828 

1.1+11+ 

1.0025 

3-61+9 

.40 

• 5203 

1.301 

.9601 

3.038 

2.2 

.22 

.3292 

1.496 

1.1466 

4.1+76 

.44 

• 5877 

1.336 

1.101 

3.589 

2.4 

.24 

.3794 

1.581 

1.2950 

5-386 

.48 

.6562 

1.367 

1.2416 

4.159 

2.6 

.26 

.4324 

1.663 

1.41+86 

6.348 

.52 

.7254 

1.395 

1.3821 

4.71*0 

3.0 

— 

— 


— 

— 

.60 

.861+1+ 

1 . 1+1+1 

1.6620 

5.913 

4.0 





— 

— 

.80 

1.209 

1.512 

2.3557 

8.821 



r — 1 

= 0.3 





^ Xt = 0.4 



1.0 

0.30 

0.3327 

1.109 

Imag. 

Imag. 

0 . 1+0 

0.1+1+02 

1.101 

Imag. 

Imag. 

1.2 

.36 

.4082 

1.134 

0.2086 

1.213 

.48 

.5370 

1.119 

Imag. 

Imag. 

1.4 

.42 

.1*867 

1.159 

.4133 

1.501 

.56 

.6358 

1.135 

0.2668 

1.420 

1.6 

.48 

.5677 

1.183 

.5854 

1.828 

.64 

.7360 

1.150 

.1+1+28 

1.677 

1.8 

-54 

.6503 

1.201+ 

• 7293 

2.184 

• 72 

.8371 

I.163 

• 5833 

1.951 

2.0 

.60 

.7341 

1.221+ 

.8636 

2.556 

.80 

.9385 

1.173 

.7100 

2.233 

2.2 

.66 

.8185 

1.240 

.9942 

2.938 

.88 

1.040 

1.182 

.8283 

2.520 

2.4 

.72 

.9032 

1.254 

1.1216 

3.326 

.96 

1.11+2 

1.189 

.9422 

2.808 

2.6 

.78 

.9881 

1.267 

1.2461+ 

3.716 

1.04 

1.21+3 

1.196 

1.0525 

3.097 

3.0 

.90 

1.158 

1.286 

1.4932 

4.1+96 

1.20 

1. 1+1+6 

1.205 

1.2508 

3.666 

3.5 

1.05 

1.369 

1.303 

1.7941+ 

5.465 

1 . 1+0 

1.699 

1.213 

1.5276 

4.391 

4.0 

1.20 

1-579 

1.315 

2.092 

6.1+23 

1.60 

1.951 

1.219 

1.7818 

5.101 


TABLE III 

VALUES OF R cr FOR A SINUSOIDALLY IDADED ARCH HAVING NONZERO AND COMPUTED FROM EQUATION (44) 

[Dashed lines indicate that there is no critical load] 

(a) m = 3 


x i 

1.0 

1.2 

1.4 

1.6 

1.8 

2.0 

2.2 

2.4 

2.6 

3.0 

3.5 

4.0 

0.01 

• 05 

.1 

.2 

• 3 
.4 

1.0000 

1.3121 

1.3079 

1.2953 

1.2547 

1.2132 

1.7615 

1.7499 

1.7161 

1.6105 

1.5011 

1.4195 

2.3637 

2.3228 

2.2493 

2.0359 

1.8285 

1.6772 

3.0881 

3.0349 

2.8926 

2.5164 

2.1840 

1.9510 

3.9955 

3.8945 

3.6492 

3.0383 

2.5562 

2.2332 

5.0877 

4.9029 

4.4765 

3.5892 

2.9384 

2.5196 

6.3812 

6.0540 

5.3864 

4.1589 

3.3260 

2.8081 

7.5720 

7.2662 

6.3485 

^.7396 

3.7159 

3.0971 

9.6960 

9.3967 

8.3545 

5.9134 

4.4964 

3.6657 

12.1039 

11.7881 

10.7127 

7.3769 

5.4650 

4.3906 

14.3783 

14.0371 

12.8865 

8.8213 

6.4228 

5.1012 


— 


(b) m = 2 


Xs \. X 1 
VM\ 

1.0 

1.2 

1.4 

1.6 

1.8 

r* 

2.0 

2.2 

2.4 

2.6 

3.0 

3.5 

4.0 

0.01 

.05 

.1 

.2 

.3 

.4 

— 

1.3121 

1.3060 

1.2887 

1.2397 

1.2017 

1.7612 

1.7422 

1.6917 

1.5616 

1.4534 

2.3476 

2.2970 

2.1790 

1.9248 

1.7328 

1.6120 

3.0837 

2.9577 

2.7243 

2.3091 

2.0236 

1.8530 

3.9775 
3 . 6889 
3.2982 
2.7019 
2.3185 
2.0910 

4.9991 

4.4435 

3.8788 

3.0964 

2.6150 

2.3306 

6.0364 

5.1904 

4.4548 

3.4892 

2.9106 

2.5705 

7.0223 
5.9191 

5.0223 
3.8788 
3.2049 
2.8102 

8.8561 

7.3218 

6.1296 

4.6475 

3.7888 

3.2878 

10.9893 

8.9975 

7.4705 

5.5905 

4.5102 

3.8814 

13.0235 
10.6163 
8.7094 
6.5179 
5.2240 
4. 4714 
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TABLE IV 


TABLE V 


VALUES OF k FROM VALUES OF R cr FOR A SINUSOIDAL ARCH 

EQUATION (62) WITH A CENTRAL CONCENTRATED LOAD 


X 1 

k 

1.0 

0.1929 x 10"3 

1.2 

• 3578 

1.4 

• 7075 

1.6 

1.408 

1.8 

2.767 

2.0 

5.486 

2.2 

10.83 

2.4 

21.70 


X-i 

^cr 

2.4 

3.089 

2.6 

3.678 

3.0 

4.716 

3-5 

5.890 

4.0 

7.000 

4.5 

8.072 

5.0 

9.122 

5-5 

10.156 

6.0 

n.179 

6.5 

12.193 

7.0 

13.201 

7.5 

14.204 

8.0 

15.206 

8.5 

16.203 

9.0 

17.195 


TABLE VI 

CRITICAL CONDITIONS FOR CENTRALLY LOADED ARCHES WITH 


Xg = 0.005X.-L AND X 3 = 0.040X.-L 


X-i 

( B 4cr 

Mcr 

Rcr 

Ncr 

C<$«r 

( a ^) cr ,max 
a P 

Max. critical 
stress 
(psi) 

(11 

1.2 

0.3713 

0.0088 

O .651 

0.075 

1.28 

3.85 

8.4 x 103 

1.5 

.6310 

.01351 

.996 

.108 

1.84 

4.50 

10.4 

2.0 

.9895 

.03263 

1.878 

.196 

2.79 

7.21 

16.3 

2.5 

1 . 4003 

.1157 

3.048 

.290 

3-57 

9.74 

21.8 

3.0 

1.9770 

•2175 

4.193 

.380 

3-73 

11.63 

25.1 

4.0 

3.0494 

.3605 

6.236 

.541 

3.78 

15.16 

30.9 

5.0 

4.039 

.4831 

8.l4o 

.693 

3-79 

18.7 

36.8 

6.0 

4.990 

• 5991 

9.986 

.842 

3.80 

22.3 

42.6 

7.0 

5.922 

.7120 

11.800 

• 990 

3.80 

25.9 

48.5 

8.0 

6.841 

.8231 

13.596 

1.136 

3.81 

29.5 

54.5 

9.0 

7.752 

• 9332 

15.380 

1.282 

3.81 

33.1 

60.3 


■''Highest outer fiber stress in arches representative of test speci- 
mens (e = 10.3 x 10 6 psi, L = 18 in., and t = 0.25 in.). 



\J1 

ON 


TABLE VII 

EFFECT OF FLEXIBILITY OF SUPPORT ON CRITICAL LOAD 


(a) Values of R cr as a function of 3 


3 \ 

1.0 

1.2 

1.4 

1.6 

1.8 

2.0 

2.2 

2.4 

2.6 

3-0 

3.5 

4.0 

1.0 

1.000 

1.312 

1.762 

2.350 

3.080 

4.000 

5.096 

6.380 

7-584 

9.708 

12.12 

14.39 

• 95 


1.288 

1.716 

2.277 

2.983 

3.850 

4.895 

6.135 

7.391 

9.564 

12.01 

14.30 

.90 

— 

1.265 

1.671 

2.204 

2.876 

3.701 

4.694 

5.873 

7.166 

9. 402 

11.88 

14.20 

.80 

— 

1.226 

1.584 

2.062 

2.664 

3.404 

4.295 

5.349 

6.584 

9.000 

11.58 

13.95 

.70 


1.200 

1.504 

1.924 

2.457 

3.155 

3-897 

4.829 

5.917 

8.439 

11.17 

13.62 

. 60 

— 

— 

1.437 

1.795 

2.256 

2.823 

3.505 

4.312 

5.255 

7.587 

10.59 

13.16 

• 50 

— 

— 

— 

1.681 

2.066 

2.544 

3.121 

3.803 

4.598 

6.565 

9.686 

12.48 


(b) Solution of equation (71) 


3 

pl)o 

1.0 

2.345 

• 95 

2. 406 

• 90 

2.471 

.80 

2.622 

■ 10 

2.803 

.60 

3.028 

• 50 

3.317 
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TABLE VIII 

VALUES OF R cr FOR A SINUSOIDAL ARCH WITH A CENTRAL 
ELASTIC SUPPORT AND q = sin ^ 


Nv M- 

0.5 

1.0 

1.5 

2.0 

3-0 

1.42 

1.6o 

1.74 

1.80 

2.00 

2.20 

2.40 

3-00 

3-50 

4.00 

2.82 

3.34 

3.87 

4.11 

5.06 

6.21 

7.43 

10.45 

17.73 

14.92 

5.10 

5.32 

6.25 

7.40 

8.43 

11.17 

13.32 

15-42 

No solution in this 
region 

7.61 

8.54 

9.38 

11.84 

13.87 

15.89 

9.28 

9.62 

10.29 

12.49 

l4.4o 

16.34 

12.46 

11.66 

12.00 

13.70 

16. 40 
17.20 
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TABLE IX 


THEORETICAL AND EXPERIMENTAL DATA 

4 


Specimen 

Length 

(in.) 

Width 

(in.) 

Thickness 

(in.) 


*2 

X 3 

Buckling 

load 

( Rcr )exp 

(Rcr) ex p 
(^cr) c i a ss 

(i) 




db) 


1 

18 

0.500 

0.249 

3-78 

0.0138 

0.136 

82.7 

5.19 

0.880 

2 


.495 

.1885 

9.12 

.0055 

.344 

85.7 

16.37 

1.004 

3 

18 

.500 

.249 

4.25 

.0097 

.146 

107.0 

6.72 

• 955 

4 

18 

.499 

.249 

3.32 

.0055 

.114 

73.7 

4.63 

.915 

5 


.501 

.249 

2.63 

.0417 

.097 

33.9 

2.11 

.653 

6 


• 493 

.249 

1.83 

.0146 

.063 

16.2 

1.04 

.671 

7 

18 

.502 

.250 

4.71 

.0842 

.159 

98.5 

6.07 

.884 

8 

18 

.502 

.251 

4.07 

.0496 

.167 

94.4 

5.72 

.976 

9 

18 

.504 

.251 

3.67 

.0666 

.164 

80.0 

4.83 

.988 

10 

18 

.505 

.250 

3.30 

.0178 

.123 

60.4 

3.70 

.781 

11 

18 

.502 

.250 

3.90 

.0264 

.126 

96.7 

5.95 

1.003 

12 

18 

.505 

.250 

5.31 

.0015 

.185 

139.8 

8.55 

.926 

13 

18 

.503 

.250 

5-07 

.0957 

.131 

115.8 

6.98 

• 925 

14 

l 7 3i 

*32 

.502 

• 374 

1.86 

.0076 

.0582 

83.3 

1.02 

.630 

15 

17 15 

17 T5 

.500 

• 375 

1.67 

.0170 

.0610 

73.0 

.886 

.703 

l6 

17— 

32 

.501 

• 374 

1.38 

.0013 

.0459 

(2) 



17 



17 15 

17 T5 

.501 

• 374 

1.265 

.0141 

.0472 

(2) 



18 

17^1 

r 32 

.502 

• 374 

2.44 

.0015 

.0850 

157.3 

1.93 

.666 

19 

17— 

32 

.499 

• 374 

2.08 

.0043 

.0707 

129.9 

1.60 

.773 

20 

17—— 
1 6 

.503 

.374 

1.34 

.0244 

.0500 

(2) 

— 

— 

21 

l 7 3i 

r 32 

.502 

• 374 

2.43 

.0112 

.0883 

176.9 

2.16 

.745 

22 

18 

.501 

.186 

6.08 

.0058 

.237 

48.2 

9.70 

.930 

23 

18 

• k 99 

.185 

6.43 

.0031 

.236 

53-6 

10.96 

.978 

24 

18 

.500 

.185 

7.23 

.0225 

.257 

62.5 

12.89 

1.031 

25 

18 

.500 

.186 

9.15 

.0007 

.311 

73-0 

14.73 

1.016 


'Material: Specimens 1 to 13 and 22 to 25, 24S-T3; specimens 14 to 21, 75S-T6. 


E = 10.3 x 10 T psi. 

p 

^Specimen did not buckle. 
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Figure 1.- Buckling mode for a high arch. 



Figure 2.- Possible buckling mode for a low arch. 



Figure 3*- Coordinate system. 


R 




Figure 4.- Relation between and R for symmetrical buckling of a 

sinusoidal arch under a sinusoidal load. 



Figure 5*- Relations between B]_, B n , and R for a sinusoidal arch 

which buckles in the -nth mode. 
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Figure 7 • - Deformation history of a sinusoidal arch. 
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Figure 8.- Change of critical load due to initial thrust H 0 . 
(ARcr) S = (Rcr) S =o * ( R cr) S= s5 S = HqL 2 /* 2 ®!. 


ON 

CO 


NACA TN 28k0 


61*. 


NACA TN 2840 




Figure 9*- Examples of low arches having non sinusoidal center lines. 


2.0 


(B 2 >cr 

X2 



(a) Ratio of (B 2) cr /^2 = C b 2) C r/ a 2* 



(b) Ratio of (B 3 ) cr /X 3 = (b 3 ) cr /a 3 . 


Figure 10.- Solution of equation (^40) for m = 2 and m = 3. 
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(a) 


m = 2 ; 


6 = 


*1 ' 


yo 

L 


sin — + ^ sin 
Xj 


mjrx 

L 


with 


m = 2 and 3- 


Figure 11.- R cr (sinusoidal) for arch forms 
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Figure 12.- R cr 


for arches under a concentrated central load. 


1.2 



Figure 13- - Variation with load of first three modes of two centrally 

loaded arches having \ 2 Ai = °*°°5 and ^Al = 0.0*10. ^ 

vo 
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Figure l4.- Effect of support flexibility on critical load. (3 = 
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Figure 15.- Critical load of a sinusoidal arch having a central elastic 

support, n = 21,3a 
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Figure 1 6 .- Theoretical and experimental results. 
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Figure 17*- Testing jig. 
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Figure 18.- Testing apparatus with specimen in place. 
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Figure 19*- Knife-edge fitting. 
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